function res = Lax_Friedrichs_onestep(u,N,mu)
u = fill_ghostcell_period(u,N,1);
res = 0.5*(u(3:N+3)+u(1:N+1)) - 0.5*mu*(u(3:N+3)-u(1:N+1));